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ABSTRACT 

The nature of the dominant component of galaxies and clusters remains unknown. 
While the astrophysics community supports the cold dark matter (CDM) paradigm 
as a clue factor in the current c osmo logical model, no direct CDM detections have 
been performed. iFaber fc Visse"rl (|2006T ) have suggested a simple method for measuring 
the dark matter equation of state that combines kinematic and gravitational lensing 
data to test the widely adopted assumption of pressureless dark matter. Following 
this formalism, we have measured the dark matter equation of state for first time 
using improved techniques. We have found that the value of the equation of state 
parameter is consistent with pressureless dark matter within the errors. Nevertheless, 
the measured value is lower than expected because typically the masses determined 
with lensing are larger than those obtained through kinematic methods. We have 
tested our techniques using simulations and we have also analyzed possible sources 
of error that could invalidate or mimic our results. In the light of this result, we can 
now suggest that the understanding of the nature of dark matter requires a complete 
general relativistic analysis. 

Key words: equation of state, gravitation, gravitational lensing: strong, gravita- 
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(cosmology:) dark matter 



1 INTRODUCTION 

Strong evidence, from a large number of independent obser- 
vations indicates that dark matter is composed by yet un- 
known weakly interacting elementary particles. Since these 
particles are required to have small random velocities at 
early times, they are called cold dark matter (CDM). Many 
solutions have been proposed to explain its presence but its 
nature remains obscure. Up to the present the hypothesis of 
pressureless dark matter remains experimentally untested 
since laborator y experiments have not yield positive results 

jBertondl201(ih. 

IFaber fc Visseil (2006) have conceived a novel approach 
to calculate the density and pressure profiles of the galactic 
fluid, with no assumptions about their specific form. Such 
test is based on General Relativity results, the weak field 
condition, and the probe particle speeds involved (photons 
and stars). In order to carry out an explicit measurement of 
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the dark matter equation of state (EoS), we have applied this 
test to galaxy clusters presenting gravitational lensing ef- 
fects. The advantage of galaxy clusters over galaxies is their 
larger dark matter concentrations and their vast spectro- 
scopic data, which allows for reliable kinematic profiles. 

2 A RELATIVISTIC EXPERIMENT 

A static spherically symmetric gravitational field is repre- 
sented by a space-time metric of the form (|Misner et al.l 
[l97oft ds 2 = -c 2 e 2 *M dt 2 +dr 2 +r 2 dQ 2 , where 
$(r) = $(r)/c 2 , and <£> is the gravitational potential. 

Resorting to the Einstein field equations with a consis- 
tent static and spherically symmetric stress-energy tensor, 
and using the mass density definition (J inp(r)r 2 = m(r)), 
the pressure profiles are: 
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Pr(r) and pt(r) refer to the radial and tangential pres- 
sure profiles, which are completely determined by the two 
functions $(r) and m(r). If these two functions are obtained 
from observations, both pressure profiles can be inferred. For 
a perfect fluid we expect p — p t = p r . 

When analyzing data, it is convenient to assume some 
simplifying hypothesis. Standard Newtonian physics is ob- 
tained in the limit of General Relativity through the follow- 



ing conditions: (i) the gravitational field is weak 
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2$ <C c , (ii) the test probe particle speeds are slow com- 
pared to the speed of light and (iii) the pressure and matter 
fluxes are small compared to the mass-energy density. While 
the first and second conditions are accomplished by galax- 
ies in clusters, photons only fulfill the first condition. The 
third condition is often applied, since it is related to the 
nature of the dominant component and to the possibility 
th at this is a pres s ureles s fluid. The novel idea introduced 
bv lFaber fc Visserl ((2006) is to avoid the assumption of the 
third condition. 

Un der condition (i) th e tt component of the Ricci ten- 
sor is |Misner et alj 119731 ) V 2 $ « -R tt , then V 2 <E> « 
^f-(c 2 p+p r + 2pt). Consequently, the function can be 
interpreted as the Newtonian gravitational potential <l?iv(r) 
if and only if the fluid is pressureless. In the kinematic 
regime conditions (i) and (ii) are fulfilled. Due to the sec- 
ond condition, the geodesic equation can be reduced to 

J2 

gpi 1 « — V<P, where r is the position vector of the galaxy 
and $(r) 7^ $jv(r). Then, the mass profile obtained from 
the kinematic analysis is defined by 
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which causes m/f(r) to differ from m(r). On the other 
hand, in the case of gravitational lensing, photons are the 
test particles and condition (ii) is not satisfied. Hence, the 
geodesic equation needs to be solved exactly to understand 
the influence of the gravitational field. By applying Fermat's 
principle and conside ring an effective refractive index (see 
iFaber fc V isscr (200^) for details), the lensing gravitational 
potential is defined as 



2$;en S (r) = $(r) + 



m(r) 



dr . 



where X7 2 &i ens (r) = 47rp iens (r). and then $>t en s(r) = 
/ m ' e "f (r) dr . This implies m iens (r) = \m K {r) + \m{r). 
This analysis gives a simple expression for the two functions 
required to calculate the density and pressure profiles: 



$(r) = 



Gm k 



m(r) = 2m iens (r) - m k . 



It is important to note that a gravitational lensing 
analysis usually assumes a Newtonian gravitational poten- 
tial, but in this general case the effective refractive index - 
physical observable of gravitational lensing- requires a more 
comprehensive definition of the gravitational potential. 

3 WEIGHING CLUSTERS OF GALAXIES 

Jeans equation offers a direct way to calculate the mass pro- 
file via kinematics: 
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where vtik (< r) is the mass enclosed within a sphere 
of radius r, p n is the 3D galaxy number density, oy is 
the 3D line-of-sight (l.o.s.) velocity dispersion and /3 is the 
anisotropy parameter /3 = 1 — {v 2 + v 2 ) /(2{v 2 )). An a lterna - 
tive approach was introduced by iDiaferio fc Gellerl l|l997h . 
who suggested the possibility of measuring cluster masses 
using only redshifts and celestial coordinates of the galaxies. 
The method they developed, called caustic technique (CT), 
allows to calculate the mass profile at radii larger than the 
virial radius, where the assumption of dynamical equilibrium 
is not valid. On a redshift space diagram (l.o.s. velocity v vs. 
projected distances R from the cluster center) cluster galax- 
ies are distributed on a characteristic trumpet shape, whose 
boundaries are called caustics. Since these caustics are re- 
lated to the l.o.s. component of the escape velocity, they 
provide a suitable measure of the cluster mass. The CT is 
a method for determining the ca ustic amplitude A ( t), an d 
then the cluster mass profile (see IDiaferio fc Gellerl l|l997h . 
and iDiaferid ((2009) for details). In order to calculate the 
kinematic mass profile we have applied the following proce- 
dure to a simul ated cluster extracte d from the Millennium 
Simulation Run (|Springel et alj2005h . and two real rich clus- 
ters of galaxies (Coma and CL0024). The galaxy systems 
were selected based on four conditions: approximate spheri- 
cal symmetry, low level of subclustering, gravitational lens- 
ing data measurements (only weak lensing reconstruction in 
the case of Coma) and an important number of measured 
redshifts of galaxies in order to accomplish our assumptions. 
Coma and CL0024 have 1119 and 271 galaxies within the 
caustics respectively. 

With the purpose of calculating the mass profile via the 
Jeans equation, we have used the first steps of the CT to re- 
m ove interlopers, an d an ad aptive kernel method (described 
in IDiaferio fc Gellerl l|l997h ) to estimate the density distri- 
bution of galaxies in the redshift diagram. In this way, we are 
able to obtain a ID profile for the l.o.s. velocity dispersion 
an and the 2D number density profile p (which are simply 
the second and first moment of the density distribution at 
each R, where R is the projected distance to the center). 
Using in this novel way the estimated density distribution 
of galaxies allows us to measure or and p at several radii, in 
order to recover the kinematic mass profile with high preci- 
sion. 

In all cases we have used a King profile to fit the num- 
ber density profile p and we have followed the procedure 
described in lDiaz et al.l |2005) to obtain the 3D number den- 
sity p n . In order to determine the velocity dispersion a r , we 
have applied the Abel inversion technique assuming f3 = 0. 
The 2D and 3D profiles are in good agreement with the 
real profiles of the simulated cluster. This indicates that the 
assumption of /3 = is quite adequate. Despite this good 
agreement, and in order to quantify the impact of /3(r) on 
the dark matter equation of state, we have solved the Jeans 
equation considering three cases: (1) /3 = and j3(r) deter- 
mined by a linear fit to (2) the data of the selected simu- 
lated cluster and (3) the data from the most massive clusters 
of Millennium, up to l/i _1 Mpc. Figure [1] shows these three 
mass profiles, together with the caustic mass profile and the 
true mass profile (NFW best fit to the simulated data) . The 
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Figure 1. Mass profiles of the simulated (left, with 1000 galaxies and left-center, with 200 galaxies panels) and real clusters (Coma, 
right-center and CL0024, right panels), for the case 1 (red dash dotted line), 2 (gray continuous line with error bars from bootstrap 
analysis), 3 (orange dash triple do tted line, only for s imulated cluster), from t he caustic technique (green dashed line) and from lensing 
analysis. Lensing p rofiles are from iKubo et al,l ll2007l) (black co ntinuous line) , iGavazzi et al,l | |2009|) (blue continuous line) in Coma and 
iKneib et al,l (|2003r l (black c ontinuous line). lUmetsu et al,l l|2010h (blue continuous line) in the CL0024. The magenta dotted li nes are the 
NFW kinematic profile from lLokas fc Mamonl 11200311 in Coma and the X-ray inferred hydrostatic equilibrium mass profile from lOta et al.l 
(2004) in CL0024. In the case of the simulated cluster, the solid cyan line is the NFW fit to the true mass profile, and the kinematic 
mass profiles were computed selecting a similar number of galaxies than those for the observed clusters, using a magnitude cut-off. 



difference between the upper panels reflects the effect of the 
number of galaxies tracing the kinematics. To mimic a lens- 
ing situation, we have projected the mass of the halos in the 
field of the simulated cluster and then we have deprojected 
the 2D density assuming that all the mass belongs to one sin- 
gle halo. The differences between the true mass profile and 
the new "fake" mass profile can be appreciated only at large 
radii and their best NFW fits are almost indistinguishable. 

In the calculation of the equation of state, we have 
combined the tangential and radial pressure, so w = (p r + 
2pt )/c 2 3p. As shown in Figure [T] there is a good agreement 
among the mass profiles of the simulated cluster (except for 
case 3), implying a null w parameter (Figure[2]). When ~ 200 
galaxies are used, the w parameter profiles differ more sig- 
nificantly. This is mainly due to uncertainties of the mass 
profile and the equation of state parameter determinations 
which originates from the low number of galaxies of this sim- 
ulated cluster. The kinematic mass profiles show also some 
variance due to the presence of inhomogeneities on the radial 
distributions of galaxies, related with the subclustering. 

In Figure [2] can also be seen that, when using the mass 
profile determined via the caustic technique, w adopts a high 
positive value in the inner regions of the clusters. This might 
happen because the CT is very effective in estimating the 
mass profile in the outs kirts but it tends to overestimate it 
within the virial radius (jSerra et al.ll20l5T ). 

4 MEASUREMENT OF THE DARK MATTER 
EOS IN COMA AND CL0024 

The results of the methods described above show a good 
agreement between the measured profile and the real mass 
profile in the case of the simulated cluster, independently 
of the number of galaxies used in the computation of the 
kinematic mass. Our method is slightly sensitive to the 
anisotropy parameter; to address this problem, we have cal- 
culated the equation of state with the kinematic mass drawn 
from cases 1,2, and 3 (explained in section [3]). The result- 
ing profiles, in Figure [2] show that the anisotropy parameter 
P has a non-negligible impact on the equation of state. As 
first result of this work we have found, in the case of the real 
clusters of galaxies, a good agreement between the kinematic 
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Figure 2. Dark matter equation of state radial profiles, corre- 
sponding to the mass profiles of fig.[T] with the same conventions. 
For clarity we split the EoS parameter determination based on the 
different mass density profiles obtained from lensing analysis. Er- 
ror bars (at 1 a level) were computed using bootstrap analysis in 
the galaxies samples used in the determination of the kinematic 
profile and in the NFW parameters of the dark matter profile 
inferred from lensing. 
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mass profiles computed using the caustic technique (heuris- 
tic recipe) and the Jeans equation inversion (rigorously cor- 
rect), but there is a noticeable difference between them and 
the mass d erived from the gravitational-lens model (Figure 
[1] see also iDiaferiol (|2009r i). However, this discrepancy al- 
lows us to measure the dark matter equation of state for 
the first time using clusters of galaxies. The resulting equa- 
tion of state of the dark matter (shown in Figure[2} behaves 
as expected when we analyze the simulated cluster, but it 
adopts an almost constant negative value for the real clusters 
(however consistent with the Strong Energy Condition) in- 
stead of the constant zero value required by CDM. We could 
attribute that to the anisotropy parameter ft. Nevertheless 
CT is not strongly dependent on j3, and the equation of state 
from the caustic profile shows also a trend to be negative. 
We have checked a number of sources of systematics such 
as the departure of sphericity and relaxation, the presence 
of deflect ing substructures in the l.o.s. (as those reported in 
Coma bv lAdami et~ail l|2009tn and the ellipticity of the halo. 
We have tested the lack of sphericity computing the mass 
profiles along three different l.o.s., showing no significant 
differences. The presence of substructures in the l.o.s. has 
been evaluated not only through the redshift distribution of 
the cluster galaxies but also deprojecting the 2D mass as- 
suming that the mass of the halos near the l.o.s. belongs to 
the same cluster (as explained in section |3J. As for the tri- 
axiality of the halo, we have modified the NFW parameters 
of the lensing mass, a ccording to the results presented in 
ICorless fc Khi3 (|2007l ), considering an extreme case of axis 
ratio Q=2.5. None o f these tests seems to explain the fea- 
tures we have shown (|Serrall2008h . although we stress that a 
combination of several of them might be responsible for this 
apparent inconsistency. 

It should be mentioned that the cluste r of galaxies 
CL00 24 experienced a merger along the l.o.s. (iCzos ke et al.l 
2002) approximately 2-3 Gyr ago. Nevertheless the good 
agr eement between t he kinematic mass profile computed and 
the lOta et alj (|2004T ) hydrostatical equilibrium mass in the 
inner region indicates that the gravitational potential has 
had t ime to relax (this was noticed first by lUmetsu et alj 
(2010)). The de nsity profiles determined by lensing m eth- 
ods in CL0024 (|Umetsu et al.ll2010l ; iKneib et al.ll2003h are 
in good agreement with each other. This is not the case of 
Coma, but it should be recalled that the Co ma density pro- 
file de rived from the weak lensing analysis of lGavazzi et al.l 
(2009) was computed in the central region (R < IMpc) us- 
ing a very deep photometry. This density profile was extrap- 
olated to outer radii in our analysis, and the w parameter 
computed using this profile shows a fair agreement w i th the 
CDM value. The weak lensing profile of lKubo et"ai] (|2007t ) 
span a much wider region using SDSS photometry, and the 
corresponding w has a constant and negative value. 

The EoS parameter recovered using this lensing analy- 
sis in Coma shows a similar behavior than those obtained in 
CL0024 (i.e. a preferred value of w ~ — |). There is, how- 
ever, a trend for t he w computed with the lensing profile 
from IKneib et al l {2003) to increase towards the external 
regions. If further measurements confirm the trend of neg- 
ative values for the dark matter equation of state w, this 
result could be interpreted in the framework of theories in- 
cluding scalar fields and the possibility of effective negative 
pressures, or alternative models of gravity. It is important to 



notice that the measured value of w in this work is consistent 
with the standard pressureless cold dark matter at la level. 
The error analysis uses 30 bootstrap resampling of the galax- 
ies in the computation of the kinematic mass profiles, and 
30 dark matter profiles resulting from the errors and degen- 
eracies of the measurements of the NFW halo parameters. 
The assumption of pressureless dark matter can be further 
tested by applying the method introduced in this work to a 
large number of lensing clusters with several galaxy mem- 
bers with measured redshifts. 
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